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1. Introduction 

The mirror Thermodynamic Bethe Ansatz (TBA) is generalization of the relativistic 
TBA [Q and offers a tool to determine the spectrum of the AdSs x S 5 superstring 
[g, Recently there has been interesting progress on its precise formulation and 
in deriving some consequences of the corresponding spectral equations |p|-[p4|. In 



parallel development the four-loop anomalous dimension of the Konishi operator 
was obtained |15[] by means of the generalized Liischer formulae [16|, [17|, [T^] 1 and 



it exhibits a perfect match with a direct field-theoretic computation ^3 |29|. More 
recently, a refined version of the generalized Liischer formula has been proposed, and 
applied to the five-loop anomalous dimension of the Konishi operator f30f , 



A (io) = A do) , 10 f _ 81C(3) 2 8103) _ 4505) 945C(7) _ 2835 } 
asympt y X 16 32 4 32 256 /' 

where g is the coupling constant (string tension) related to the 't Hooft coupling 
A through A = 47r 2 g 2 . This result has been further generalized to all twist two 
operators by also invoking the reciprocity principle [[31 . Although no field-theoretical 
computation has been done at five loops, those anomalous dimension being continued 
to the negative values of spin enjoy a quite non-trivial agreement with the constraints 
imposed by the BFKL equation p2|-H35 . 



Having in mind all these developments, it is time to ask whether the proposed 
TBA spectral equations are in accord with the perturbative findings based on the 
generalization of Liischer's approach. As for Konishi at four loops, incorporation 
of the corresponding result in the TBA framework does not pose any real difficulty 
because the leading finite-size correction has been built-in when the excited-state 
TBA is formulated; all the contribution to the anomalous dimension (the energy 
of the corresponding string state) comes from the main Ig-functions taken at their 



1 See ffjl-[pO| for other interesting applications of Liischer's approach. 
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asymptotic values ||, the latter being given by the generalized Liischer formulae 
constructed through the infinite- volume scattering data |]30| . Also, rapidities Uk, 



k = 1, . . . , N, of the excited string particles forming an iV-particle state under con- 
sideration are determined by the asymptotic Bethe Ansatz (ABA) equations (also 
known as the Bethe- Yang equations). 

The situations change, however, for the five-loop case. As was argued in [15|, |30|| , 



to find the anomalous dimension of Konishi at five loops, one has to compute the 
correction to the ABA, though still the asymptotic Yg-functions can be used. The 
exact Bethe equation determines the shift 5uk of particle rapidities Uk = u° k + g 8 5uk 
from their asymptotic values u° k as 



M 



5kBk(u k ) 



The first term here is a variation of the ABA equations and the term $^ 8 ^ is the 
leading finite-size correction to the ABA of order g 8 . The formula above should be 

(2) (2) 

evaluated at Uk — u k , where u k are the one-loop values of the particle rapidities. 
As soon as 5u k are determined, the five- loop correction to the dimension (energy of 
the corresponding string state) follows from expanding the exact dispersion relation 
up to the five-loop order, i.e. up to g 10 . 

It turns out that the correction $( 8 ) derived through analytic continuation of the 
mirror TBA equation looks rather different from the one derived through Liischer- 
type perturbative arguments p0| . The difference occurs due to the fact that the 



exact Bethe equations involve auxiliary Y-functions which must satisfy a coupled 
system of the TBA equations. Starting from five loops, the auxiliary Y-functions 
start to contribute non-trivially to the modification of the ABA. 

In this note we find a strong evidence that the results obtained for the Konishi 
operator from both the TBA and Liischer approaches are in an excellent agreement 
at the five-loop level. Fortunately, as discussed above, what we have to do is to show 
only that the correction to the ABA derived from the TBA agrees with found 
from Luscher's approach j3(J. This will be done by expanding the TBA equations 
around the asymptotic solution || and by linearizing the exact Bethe equations 
around the ABA. As we will see, the exact Bethe equations at order g 8 involve 
a leading correction to a single auxiliary function Yx\ vw , which we will determine 
from the linearized TBA equations numerically. Then, using this result for Y\\ vw we 
evaluate numerically the correction to the exact Bethe equations and find that it 
agrees with $( 8 ) of |HJ with a sufficiently high precision. It would be important to 
support this numerical agreement by an analytic proof. 

The paper is organized as follows. In the next section we present the linearization 
of the simplified TBA equations. We use them in section 3 to compute numerically 
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the correction to the exact Bethe equations finding an agreement with the corre- 
sponding result in [jBOfl . 



2. Linearizing the TBA system 

In the light-cone gauge the string vacuum corresponds to the gauge theory operator 
trZ J , where Z is one of the three complex scalars of the Af = 4 super Yang-Mills 
theory. In this work we will be interested in the excited states from the s((2)-sector, 
the latter comprises the composite operators of the type tr D N Z J , where D is a light- 
cone derivative. For J = 2 these are operators of twist two spin N. In particular, 
the operator with J = 2 and N = 2 is the sl(2) descendent of the Konishi operator. 

The Thermodynamic Bethe Ansatz approach for the AdSs x S 5 mirror model 
leads to the following expression for the energy of the corresponding iV-particle string 

states 

E = J + J2S( Pl )--J2 d«|-l g(l + y Q ). (2.1) 

i=l Q=i J -°° 

Here J is the angular momentum carried by the string rotating around the equator 
of S 5 . The integration runs over a real rapidity line of the mirror theory, and Yq 
are momenta and Y-functions of the mirror Q-particles. The asymptotic energies of 
string theory particles with momenta pi are fixed by the dispersion relation [R6|, E7\ 



S{p) = A /l + Vsin 2 ^. (2.2) 



The function Yq is exponentially small at large J. Therefore the last term of ( |2.1| ) can 
be regarded as a finite-size correction to the asymptotic spectrum of string energies. 

For a fixed J and small g, the Yq functions become small again. The finite-size 
corrections provide so-called wrapping corrections to the energy or the anomalous 
dimension at weak coupling. In particular, for the Konishi operator the finite-size 
effects make their appearance starting from g 8 that corresponds to the fourth loop 



order of perturbation theory [15 |. 



Recall that the large J asymptotic solution of the excited-state TBA equations 
can correctly reproduce the leading finite-size corrections to the energy. The TBA 
equations are formulated in terms of the following Y-functions: Y^-functions associ- 
ated with Q-particle bound states, auxiliary functions Yq\ vw for Q|t>w-strings, Yq\ w 
for Q|u>-strings, and Y± for y-t-particles [|} |j. 

In what follows it is convenient to use the simplified TBA equations derived in 
[| in]. To determine the leading finite-size correction to the asymptotic form 
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of Y-functions, we linearize the simplified TBA equations by introducing for any 
Y-function the following representation 



Y = Y°(l + #0, (2.3) 

where Y° is the corresponding asymptotic expression and *3f is treated as the per- 
turbation. The linearized TBA equations will then take the following form 

• M|w-strings: M > 1 , %\ w = 

&m\w = (Am-i\w&m-i\w + A M+ i\ w <& M+l \ w ) * s + 6 mi \ i_y° ~ i - y° ) * S ' ^ 2 ' 4 ' ) 

Y° 

where we have introduced the concise notation A M \ W = , - h ' £ - . 

• iW>w-strings: M > 1 , %\ vw = 

&M\vw — (Am-1\vw^M-1\vw + Am+1\vw&M+1\wo) * s ~~ ^M+l * s (2-5) 

+ 6m1 ( 1 _ A. ~ 1 _^ J* g ' 

^ 1 Y° y° 

Y° 

Here we also defined A M \ VW = . Any asymptotic Y-functions except for 

1 1+r M\vw 

Yq have the magnitude of order 1, so A M \ W , A M \ VW are not small. Note that the Y- 
function provides the leading large J correction to the asymptotic TBA equations 
and, for this reason, it enters the last equation as an inhomogeneous term. 

• ^/-particles 

& + -&L = Y$*K Qy , (2.6) 
& + + &L = 2(A 1]VW & 1]VW - A 1{w & 1]w ) *s — Yq*s + 2Yq * K® 1 * s . (2.7) 

We will not need the equations for Q-particles, because the asymptotic solution 
for Yq is already known from the generalized Liischer's formulas. The equations for 
Q-particles suggest that the next corrections to Yq start from g 16 . 

The linearized TBA equations above define the leading, exponential in J cor- 
rection to the asymptotic Y-functions. However, by further fixing J = 2 we would 
like to view W as another source of perturbative correction in g to the asymptotic 
Y-functions. Therefore it makes sense to expand further the asymptotic Y-functions 
in the linearized TBA equations in powers of g having in mind that expansion of Ym 
starts from g 8 . 

The rapidities Uk of the Konishi operator are expanded at small g as u% = — w 2 — 
1 + 2 g 2 — + 0(g 6 ) using the notation of [14]. It is convenient to rescale 



l 

V3g 



them as u& — > Uk/g and do the same u — >■ u/g with the argument u of the Y^(u) 



- 4 - 



Ym\vw( u ) anc ^ ^m\w( u )- ^ n ^ ne other hand, the functions Y±(u) are supported on the 
segment u G [—2, 2] and as for them, the variable u will be kept unrescaled. 

First, let us consider the equation for y-particles. Performing the change of 
variables u — > u/g in the convolution term Yq * s and expanding in g, we see that 

Y°*s- r dU Y°( U ) 9 r du Y % {U) (28) 

Q loo 9 rQ l^4coshf(H-|)- J /_ oo d %coshf( M -,)' (2 ' 8) 

where the rescaled Y-function Yq(u) reads as 

Y°( u ) = a 8 64Q 2 [-l + Q 2 + n 2 - W 2 ] 2 

1 



10 A 



+ 0(g 

and we denote by w = U\ = —u 2 = 4= + 0(g 2 ) the rescaled rapidity of the two- 



[(Q - l) 2 + (« + u;) 2 ][(Q + l) 2 + (u + wf] 
i 



particle state corresponding to the Konishi operator. Thus, Yq * s starts at order g 8 
in the weak-coupling expansion. Analogously, one can establish that in eq. (|2.6| ) at 
g —7- one has Yq * Kq v ~ g 9 , while in eq.( [2.7|) one finds Yq * * s ~ g 8 . As a 
consequence, eqs. fl2.6Q and Q2.7|) imply that the first non-trivial corrections ( 3f± start 
at order g 8 , although for the difference one has — ~ g 9 . 

In fact, there is a general pattern behind inherited from the asymptotic solution 
for Y± — the sum + admits an expansion in even powers of g, while the 
difference — *3C in odd. Another interesting property is that &+(g) = *3C(—g). 
Thus, in the perturbative expansion, the coefficients of ( 3f± in front of even powers 
of g coincide. 

Now we turn our attention to eq. Q2.5l) . The notation £ signifies that integration 
in the corresponding convolution term is over the segment [—2,2]. Due to the fact 
that we do not rescale the integration variable in the integrals involving the 
expansion of s in this convolution term starts from g. Moreover, one finds that at 
the leading order in g: 

w 2 - 3 

Y ° = Y - = ^T-v (2 - 10) 

which is a w-independent quantity. Hence, 

( W - ^ + -)*s 9 =\{l + w 2 ){^-& + )*s~g w 



y, _ j_ 1 _JL J /" U 4 V 1 ;v + > y ' 

and, therefore this term produces higher than the leading g 8 contribution into the 
equation for &~m\vw 

Analogous consideration can be applied to eq. ([2.4|) which shows that the cor- 
rection &m\w starts from the order g 10 . Thus, at order g 8 we have the following 
equations to be considered 
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Figure 1: In the upper picture the profiles of the functions &i\ vw /g 8 (blue) and ^^w/^ 
(purple) are depicted. They are obtained by solving numerically eq.(2.13). The lower 
picture contains the profiles of &m\vw/9 8 f° r M = 2,3,4,5. The absolute value of &m\vw 
decreases as M increases. 

• iW>w-strings: M > 1 , %\ vw = 

&M\vw — {■^-M-\\vw ( 3^M-\\vw + Am+1\vw&M+1\vw) * s ~ ^M+l * s - (2-H) 



• ^/-particles = 

2& + = 2A l[vw & llvw *s-Y$*s + 2Y$* Kg 1 * s , (2.12) 

where all the TBA kernels are taken at their leading order. In fact, as we will see 
below, the correction to the asymptotic Bethe-Yang equations involves at order g 8 
the functions ^m\vw alone, and therefore, the only equation we have to solve is (|2.11| ). 
Eq. fl2.1ip can be equivalently written as 

%\v W * Sl QM = -Y° i+l * 8 . (2.13) 

Here we denned a kernel 



Oqm(m, v) = 6 QM S(u - v) (2.14) 
- 8q,m-iAm-i\vw{u) s i u -v)- S Q} m+iA m +i\vw(u) s(u - v) , 



where 



*M\vw 



(u) 



M(M + 2) 
(M + l) 2 



x 1 + 



1 + w 2 



1 + w 2 



w 2 v? + p 2 + 



y/M(M + 2) 

with p| = (1 ± v/M(M + 2) - w 2 ) 2 and s(it) 



y^M(M + 2) 



w 2 u 2 + 



(2.15) 
+ 0(5 2 ), 



4 cosh 2^ • 



The formal solution to 



eg. ( |2.13| ) can be given with the help of the inverse kernel Q 1 



(2.16) 



The kernel Q can be inverted by the power series expansion and computed numeri- 
cally. This gives rise to a numerical determination of the correction < 3f\\ vm . Alterna- 



tively, eq. Q2.13| ) can be solved by iterations, and it is what we have used. 
We studied the following inhomogeneous linear equation 



M\vw 



(2.17) 



which can be solved by numerical iterations starting with 



M\vw 



for each n 



2, 3, . . .. Since all equations are linear, the solution of the original equation ( p. 13 ) is 

m\vw T ne solution is unambiguously well-defined 



given by the sum 1V M \vw = Yln=2 ^» 



provided that the homogeneous equation has only the trivial solution ^m\vw 
The profiles of the first few &m\vw are presented in Figure 1. 



0. 



3. Correction to the asymptotic Bethe equations 

The exact Bethe equations determine the exact (non-asymptotic) positions of the 
Bethe roots u k . As was shown in [TA\, for Konishi-like states and below the first 
critical point the exact Bethe equations 2 admit the following representation 



N 



n(2n k + 1) = Lp k + i ^ log (uj, u k ) + Im U 



(3.1) 



2 J2 Ml + Y M ) * (lmKZ u - a*Im Jf^ 1 ' 1 -) 



M=l 



+ 2\og(l + Y llvw )*(lm§-is*K yU ) 



2 log *s *hnK% x -zlog (l - i-) (l - ^Ui^ 1 * . 



2 According to 38, 3i|, the exact Bethe equations are Yi.*(itfc) = — 1, where Yi* is the function Y\ 
analytically continued to the real rapidity line of the string region. The exact Bethe equations which 
follow from the canonical TBA equation for Y\ have been discussed in p2| . We exploit here the 
exact Bethe equations in a different representation which is derived from the hybrid TBA equation 
for Y\ discussed in (H. 
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The kernels Kf n , and Kf^* with one leg on the real rapidity line of the 



mirror theory and the other one on the real line of string theory are given in [fL4 



The quantity U appearing in the first line of the exact Bethe equation is defined as 



U(v) 



/+oo " 
at log n%-#-%) 

3=1 

N 

-2 M«j 



Finally, the kernel s is defined as 



2 2 n 

v ) — 

9 x* 



(3.2) 



s(u, v) = s{u — v h z'O) . 



(3.3) 



The second variable of any kernel entering eq.( [3.1|) and the function U(v) should be 
evaluated at v = u k - Only the imaginary part of the exact Bethe equations is written 
down in eq. ( |3.1|) , because using the explicit form of the kernels one can recognize 
that the real part vanishes, as it should be for the equation to have a real solution. 
Note that eq.( |3.1|) does not involve 1m|«t functions at all, and out of Y M \ VW only Yi\ vw 
is present. 

The exact Bethe equations can be written as 



N 



ir{2n k + 1) = Jpk + i^hg Sj^J (uj, u k ) + TZ k 



(3.4) 



Now we would like to apply the same strategy as the one in the TBA equations and 
find the leading weak coupling correction to the asymptotic value 7^ ympt = o. The 
condition 7^ s y mpt — o must be satisfied for any g to guarantee the compatibility of 
the exact Bethe equations with the asymptotic Bethe ansatz. Taking into account 
that for the Konishi state J = 2 and L = 4, we have checked numerically that this 
is indeed the case, for arbitrary values of w = u± — —u 2 inside the region where only 
Ym vw has two zeroes in the strip \\mv\ < l/g. Thanks to this identity, we can remove 
the term 5wd w ll from the corrections to the asymptotic Bethe ansatz. 

Next, one can see that the leading weak-coupling correction, which we denote by 
TZ k starts at g 8 . It turns out that the last two lines in eq. Q3.1p contain the convolution 
* which, just as in the case of the TBA equations, will not contribute to the leading 
order of perturbation theory. Thus, at leading order we have the following shift of 
the asymptotic Bethe- Yang equations 



-2 J2 Y M * Im 4i, + 2A 1{vw & 1{vw *Ims + 2j2 Y^ +1 *s*Im K^* x . (3.5) 



M=l 



M=l 
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Here all the kernels are evaluated at g = which provides the leading order contri- 
bution. In particular, we have 



lmK ll*x = -77" 

Z7T 



U — V 



1 



4+ (u - f) 2 



p.v.- 



u — V 



(3.6) 



The convolution of the kernels in the last term of eq. (|3.5|) can be computed in terms 
of Digamma functions; for all M > 1 we find 



s * Im K, 



MU 

vwx 



1 

$717 1 



— 1 



4(« - v) 



M 2 + (u- v) 2 
+ V>(f + {(u -v))-1>{f- \{u - v)) 

-^(^+i(«- v ))+^(^-3(«-«)) 



(3.7) 



Further, one needs Ims 



Finally, the leading contribution of the 



P- V '4smhf (u-v)- 

imaginary part of the dressing kernel in the mirror-string region is given by 



u — v 



2vr[(M + l) 2 + (u - v) 



(3. 



The $ and ^ functions which appear in the definition of the dressing phase kernel 
[ PH] , do not contribute at this order, because they can be written as integrals over 
the interval u 6 [—2(7, 2g] after the rescaling of rapidities. Combining now everything 
together, we read off the necessary correction 3 



00 

E 

M=l 
+ 



dU Ym(U) —- — : — 

AA ' 7T (M + l) 2 + (U - U k ) 2 



duA llvw (u)^i\ vw (u) 



M=l 



-00 
00 



2 sinh I (« — u k ) 
. 4(u-u fc ) 



(3.9) 



M 2 + (m - M fc ) 2 



Af+2 



Mfc)) 



In |!0| the leading correction to the asymptotic Bethe-Yang equations was de- 
noted by $( 8 \ Comparison of our definitions with that of shows that the agree- 
ment between the Liischer and the TBA approaches relies on the fulfillment of the 
following equality 



5K k + & 8 \u k ) = 0. 



(3.10) 



3 In the second line of the equation below the function l/sinh(w — u k ) exhibits singularity at 
u = Uk- We can however omit the principal value prescription for the integral because A^ vw 
vanishes at u — u k making the integrand regular everywhere. 
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w 



-j= = —u 2 is the one-loop rapidity for the two-particle state corre- 
sponding to the Konishi operator. 



where u\ 
nding 
According to 

$( 8 ) 



the quantity $ < - 8 - ) was found to be 



u k ) = / duY Mi u ) x 

M=l 

1 r u — Uk 



(3.H) 



x 



7T 



U — Uk 



+ 



u k 



(M+ l) 2 + {u-Ukf (M - l) 2 + (u - u k ) 2 -l + M 2 + u 2 -u 2 k 



Apparently, formulae eqs.( |3.9| ) and ( |3.11| ) look rather different, only the first line in 
eq.( |3.9| ), which is due to the contribution of the dressing kernel, cancels in the sum 
with the first term in eq.( |3.11| ). 



Plugging in eq. (|3"D|) the representation ( j2.16[ ) specified for &i\ vw , we see that the 
last two lines in ( p.9[ ) never involve Y°. The last two terms in eq. (|3.11|) exhibit the 
same phenomenon: for M = 1 they cancel each other. This suggests to compare 
individual terms in eqs.( [3~9l , [3.11|) associated with a contribution of a given Ym- The 
case of M — 1 trivially matches, and the term with &i\ vw should be understood 
as defined in ( [2.17|) . Performing numerical evaluation of the corresponding 

quantities, we summarized the results in Table 1 for M — 2, ... ,6. 



Table 1. Comparison of the numerical results 



M 


TBA 


Liischer 


TBA-Luscher 


TBA-Liischer 
Liischer 


2 


-0.0108303 


-0.0108304 


5.21829 • KT 8 


4.81819 • lO" 6 


3 


-0.000118621 


-0.000118665 


4.34822 ■ 10~ 8 


0.000366429 


4 


-4.63638 • 10~ 6 


-4.6417 ■ 10~ 6 


5.32116 • KT 9 


0.00114638 


5 


-3.78671 • 10~ 7 


-3.79407- 10- 7 


7.35899 • 10- 10 


0.0019396 


6 


-4.89345 • 10~ 8 


-4.94077 • 10~ 8 


4.73166 • lO" 10 


0.00957676 



The table describes a numerical comparison between blZx and — ^^(iti). In the 
column "TBA" the results for the contributions to 51Z\ produced by individual 
Ym's are presented for M = 2, . . . , 6. Analogously, the column "Liischer" represents 
the same quantities but for — <3?( 8 )(ui). The contributions of the first line in fl3.9|) 



and the first term in (3.11) are not included in the table. The other two columns 
give the absolute and relative estimates for the difference between the TBA and 
Liischer approaches. 
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Albeit looking different, the individual Y^-contributions show a good agreement 
between 511 k and — & 8 >. Increasing M the precision tends to decrease, but this is 
because the contributions corresponding to higher M become rather small and their 
evaluation requires more time. 

We conclude this paper with an observation to be understood in future. On the 
one hand, the expression (|3.11| ) can be concisely written as 

^\ U1 ) = -^\ U2 ) = -J2 / (3.i2) 

M=1 J-oo 47r 

by using ( p.9|) . Here the level-matching condition is imposed before the differentia- 
tion. The right hand side involves the derivative of the transfer matrix, and looks 



similar to the known result for the 0(4) model [41]. On the other hand, the gen- 



eralized Liischer formula proposed in |30] involves the trace of the derivative of the 



S-matrix with respect to the mirror rapidity u. The two quantities coincide for the 
present problem of Konishi at five loops, but must be distinguished in the general 
situation. 4 
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